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Resonances and final state interactions (FSI) play a role in the formation of CP violation (CPV) 
constrained by CPT invariance. We provide a general formulation of CPV including resonances and 
FSI starting from the CPT constraint. Our discussion is elaborated within a simple B decay model 
with the p and /o(980) resonances plus a non resonant background including the tttt —>■ KK coupled 
amplitude. We consider few illustrative examples to show the interference patterns appearing in 
the CP asymmetry, namely, that from the p resonance plus a non-resonant amplitude, and that 
from the interference of the p and /o(980) resonances. We perform the fit of the CP asymmetry 
for the charmless three-body decay channel B^ —>■ and obtain as outcome the B^ —>■ 

tv^K+K- for TTTT channel asymmetry in the mass region below 1.6 GeV in fair agreement with the 
new data LHCb data. Analogously, we also describe the CP asymmetry of the B^ —>■ 
decay, with that from the B^ —>■ K~ channel obtained as output. As in the previous case, 

we also found agreement with LHCb experimental data. 

PACS numbers: 13.25.Hw, 11.30.Er, 11.80.Gw, 12.15.Hh 


I. INTRODUCTION 

The exact CPT invariance implies the identity of the 
lifetime of a particle and its charge conjugate. Therefore, 
when the individual partial decay widths of CP conjugate 
channels are different, due to the CP violation (CPV), the 
other channels must have equal amount of CP violation, 
with opposite sign, such that the total width of the par¬ 
ticle and antiparticle are equal. Final state interaction 
must be the responsible to distribute the CP asymmetry 
among the different conjugate decay channels such as the 
sum of partial widths provides identical total width for 
the particle and its anti-particle. 

The common belief says that CPT invariance is not 
a practical constraint to be taken into account when 
computing CP violation in charmless B decays. It is 
naively expected that CP violation is distributed among 
several different coupled final state channels with two, 
three, four and more hadrons [1]. Therefore, the CPT 
constraint should be realized, in principle, by summing 
up over many hadronic decay channels. However, since 
hadronic many-body rescattering effects are far from be¬ 
ing understood, it is evident that this phenomenological 
hypothesis should be better understood before consider 
its important consequences for the description of the CP 
violation phenomena. 

In an opposite direction to the common believe, we 
shown in detail in a previous paper [2] , that the coupling 
between several final states must be suppressed at least 
in some charmless three body B decays. A quick ob¬ 
servation to the event distribution of these decays shows 
that they are placed basically at the boundary of the 
Dalitz plot, in a region dominated by light two-body res¬ 
onances. This dominance can be directly associated with 
the non-perturbative regime at the B mass energy [3] . In 


addition, by considering the observations made by early 
experiments [4, 5] that the low mass tttt and ATtt region 
is dominated by the elastic regime, the number of cou¬ 
pled channel expected for charmless three-body B decays 
should be very restricted. There is a significant excep¬ 
tion to this general picture for the two-body mass distri¬ 
bution: the S-wave tttt —>■ KK re-scattering process be¬ 
tween 1 to 1.6 GeV [6, 7], that was also discussed in detail 
in our previous paper. So, at least for this class of de¬ 
cays, the existence of few possibilities of coupled hadronic 
channels demands the explicit consideration of the im¬ 
portant constraint from CPT invariance in the study of 
CP violation. Unfortunately little is known, either theo¬ 
retically or experimentally, about the event distribution 
for others coupled channels with four or more hadrons 
in the final state. However, by looking at the profusion 
of observed low mass vector-vector charmless B decays, 
it seems that they are also dominated by low mass reso¬ 
nances [8]. 

The inclusive CP asymmetry recently released by the 
LHCb collaboration for four charged B charmless three- 
body decays, gives a hint about the correlation among the 
asymmetries in the coupled decay channels [9] . Actually, 
the experimental results for strangeness A = 1 decays 
Acp{B^ AT^tt+tt-) = -t0.025± 0.004±0.004±0.007 
and Acp{B^ K^K+K-) = -0.036±0.004±0.002± 
0.007 multiplied by the observed number of events [9] and 
corrected by the branching fractions (see Ref. [8]), which 
represent approximately the efficiency not reported in 
the paper, show that the number of events involved in 
the CP asymmetry in one decay channel is approxi¬ 
mately the same in the other channel but with oppo¬ 
site sign. The same happens, within the experimental 
errors, to the coupled decay channels with strangeness 
A = 0. Taking into account the experimental values 
Acp{B^ 7r±7r+7r-) = ±0.058 ± 0.008 ± 0.009 ± 0.007 
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and Acp{B^ n^K+R-) = -0.123 ± 0.017 ± 0.012 ± 
0.007, the number of events observed in each channel, 
corrected in the same way as before by the branching 
fractions, show the same correlation as above, namely, a 
similar number of events involved in the CP asymmetry 
but with opposite signs. 

Beyond the inclusive CPV, three-body decays allow to 
observe this asymmetry in the Dalitz phase-space. In 
principle, we can expect CP violation in charmless three- 
body B decays coming from three different types of in¬ 
terferences involving weak and strong phases: (i) direct 
CPV due to the interference of the tree and penguin am¬ 
plitudes in the same intermediate state. This first kind 
of CPV is based in the BSS model [10]. (ii) CP asym¬ 
metry produced through the interference between two 
different final states with different weak phases coupled 
by the final state interaction. This second case is con¬ 
strained by the CPT invariance and it was already well 
described in our previous work [2] . (iii) CPV from the in¬ 
terference between two neighbor resonances in the Dalitz 
plot, which share the same-phase space region. This third 
CPV kind has been pointed out in previous theoretical 
papers [11-16]. 

The present study addresses the three types of dy¬ 
namics related with CPV in the Dalitz plot discussed 
above and, more specifically, we look for signatures in¬ 
volving each one of these contributions in the charm¬ 
less decay channels B^ —>■ 7r=*^7r+7r“, B^ tt^K+R-, 
B^ —>■ iir^7r+7r“ and B^ —>■ R^R'^R~. To perform this 
study, we use the recently published LHCb paper [9], in 
particular, the tt+tt” and R'^R~ mass distribution from 
the difference between B~ and B+. These differences 
other than give directly the CP violation distribution 
in the phase-space, minimize substantially problems of 
acceptance and background not available in the LHCb 
paper. 

Our work develops a CP asymmetry formula including 
resonances and FSI based on the CPT constraint. The 
CP asymmetry is derived in lowest order in the strong 
interaction and decomposed in angular momentum. The 
explicit expression for the CP asymmetry is found for the 
decays channels B^ —>■ and B^ —>■ R^tt^'k~. 

Our discussion is exemplified by considering, within the 
isobar model, the interference patterns for the CP asym¬ 
metry from the p and /o(980) resonances plus a non 
resonant background including the contribution of the 
TTTT ^ RR coupled channel amplitude. We use this for¬ 
mula to fit the recent B^ —>■ 7r=*=7r+7r“ LHCb data [9] and 
obtain as outcome the B^ —>■ tt^R'^R~ in qualitative 
agreement with data for the mass region below 1.6 GeV. 
We apply the same procedure to the B^ —iF^7r+7r“ 
and B^ —>■ R^R'^R~ decays. We also take into account 
the LHCb separation between positive and negative cos 9 
values, that supplies new details which reveal the role of 
the vector meson resonance in building the CP asymme¬ 
try interference pattern. 

The work is organized as follows. A brief review of 
the CPT constraint for deriving the CP asymmetry with 


FSI and our notation are provided in Sec. H. The CP 
asymmetry formula in leading order of the strong inter¬ 
action and angular momentum decomposition is given 
in Sec. HI. The resonances are introduced in the CP 
asymmetry expression in Sec. IV. The formula for the 
interfering resonant and non-resonant amplitudes in the 
asymmetry is developed in Sec. V. The analysis of the 
various CP asymmetry terms is given in Sec. VI. The 
results from the fitting procedure to the asymmetry in 
the B^ —^ TT^TT+TT are shown in Sec. VH, where we 
found that the CPT violating terms can be disregarded. 
We also show results for B^ —)► 'k^R~^R~ decay with 
no new parameters. An analogous study is performed 
at Sec. VHI for B^ —>■ R^tt~^tt~ and B^ R^R~^R~ 
decays. The concluding remarks are provided in Sec. IX, 
where we also discuss a more general formula to the CP 
asymmetry, where the FSI terms interfere with the p 
and /o(980) resonances. In the Appendices, we show the 
angle integrated asymmetry as well as the derivation of 
kinematical factors. 


II. CPT INVARIANCE IN A WEAK DECAY 

We follow closely Refs. [2, 17, 18] to introduce our no¬ 
tation and the CPT constraint in B meson decays. We 
start by reminding that the weak and strong Hamiltoni¬ 
ans conserve CPT, namely, 

{crry^H^ (h,) cvt = h^ (h,) . (i) 

The hadron weak decay amplitude is computed from the 
matrix element {Xout\Hu,\h), where the distorted state 
Xout has the effect of the hadronic strong force due to 
FSI. The requirement of CPT invariance on the decay 
amplitude is 

{Xout\Hw\h) = XhXx{Xin\Hw\h)* , (2) 

where we used that the hadron state \h) transforms under 
CPT as CVT\h) = x(/ij, where h is the charge conjugate 
state and x ^ phase. 

The completeness relation of the strongly interacting 
states, eigenstates of Hg, and the hermiticity of im¬ 
plies that 

{Xout\Hw\h) = X/iXa x(AQ„f[i7m|h) , (3) 

a' 

where 

%',A = = 5v.a (4) 

defines the matrix elements of the S-matrix. The sum of 
the partial widths of the hadron decay channels and the 
correspondent sum for the charge conjugate are identi¬ 
cal, i.e. the mean life of the particle and its conjugate 
are equal, which is a consequence of Eq. (3) and the her¬ 
miticity of Hyj (see Ref. [2]): 

Xk-' out |F7ii;jh)[ —1 (Aoui[Hujjh)[ . 

A 


( 5 ) 
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The CP-violating phase enters linearly at lowest order 
in the hadron decay amplitude as suggested by the BSS 
mechanism [10]. In general, the decay amplitude can be 
written as 

A^=Ax + ( 6 ) 

where A\ and B\ are complex amplitudes invariant under 
CP, containing the strongly interacting final-state chan¬ 
nel, i.e., 

A- ={\out\H^\h), (7) 

and 

A+= {Xout\H^\h)- (8) 

The only change due to the CP transformation is the 
sign multiplying the weak phase 7 . The CPT condition 
depicted by Eq. (5) gives 

Er(^A) = Er(^^)’ ( 9 ) 

A A 

where the subindex labels the final state channels, 
summed up in the kinematically allowed phase-space. 
The decay amplitude written in terms of the CPT con¬ 
straint (3), and considering the CP violating amplitudes 
for the hadron and its charge conjugate is 

Aa + e^^-^SA = XhXx E (10) 

A' 

Note that this equation imposes a relation between A\ 
or B\ with their respective complex conjugates. 

III. CP ASYMMETRY AND ESI IN LEADING 
ORDER 

The S-matrix is given in terms of the scattering ampli¬ 
tude tA',A, namely, 

•^AyA = ^A',A + * ^A'.A, (11) 

which turns Eq. (10) into 

A;, + = XhXx {Ax + e^^'^Bx)* 

+ iXhXx^tx',x i^x' + e^’^'^By)* ■{^‘2) 

A' 

We now consider effect of the ESI in leading order in 
the decay amplitudes Aa and Bx- For that purpose, let us 
introduce decay amplitudes computed without the effect 
of the final state interaction, 


{Xo\Hyj\h) = Aqa + e “"’'Boa, 

(13) 

{Xo\Hw\h) = Aqa + e“’'“"’'i?oA, 

(14) 


where |Ao) and jAg) are the mesonic noninteracting 
charge conjugate states, i.e., without the distortion of 


the strong hadronic interaction. The terms Aqa and Bqx 
can be in principle associated with the tree and penguin 
amplitudes in the BSS model [10]. 

The leading order (LO) effect due to ESI in the decay 
amplitude is obtained by substituting Aa —>■ Aqa and 
Bx —?> Bqx in the right-hand side of Eq. (12). Consider¬ 
ing, the assumption of the CPT invariance of the weak 
Hamiltonian one has also that: 

(XolHjh) = XhXx(Xo\Hjh)*, (15) 

which implies in the following relation for the partonic 
amplitudes 

^OA = XhXxAo^ (16) 

and 

Box = XhXxBlx, (17) 

and therefore up to the leading order in tA'.A, Eq. (12) 

reduces to 

A^o = ^OA + e^^-^HoA 

+ z^fA'.A(AoA'+e±*^Hov), ( 18 ) 

A' 

which is equivalent to the one provided in Refs. [19, 20]. 

In lowest order the scattering amplitude is restricted 
only to two-body terms, therefore it is useful to decom¬ 
pose the source decay amplitudes Aqa, Hqa and the Ia'.a 
matrix in angular momentum states J of each pair in 
the outgoing channel. The angular momentum decom¬ 
position will allow to identify and introduce the meson 
resonances, like the p(770) or /o(980) in the tttt channel. 
This leading order decay amplitude (18) decomposed in 
J is written as, 

^LO=EKA + e±*'^ilj(A) 

J 

+ * E ^a'.a (^oa' + e^*"’'71oA') ) (19) 

X'.J 

where one should note that A (A') refers to two-body 
channels in the final hadronic state including all other 
dependences on the quantum numbers and in the energy- 
momentum of the spectator hadrons. By using this ex¬ 
pression, the CP asymmetry can be written as 


M'x=T [h ^ X) - T[h ^ X) 

= 4(sin7) ^Im{(Ho"';,)* A^^ 
J J' 



where the A' represents each state coupled by the strong 
interaction to the decay channel A. We just remind 
that the second and third terms in the right-hand side 
of Eq. (20) can be associated to the “compound” CP 
asymmetry [21]. These two terms cancel each other when 
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summed in all channels A and integrated over the phase- 
space, that leads to the CPT condition expressed by 
Eq. (9), once the source term in (20) satisfies 


^Im 

A J 




= 0 , 


( 21 ) 


which is a consequence of the CPT constraint at the mi¬ 
croscopic level, e.g., as expressed by the tree and penguin 
amplitudes in the BSS model, that should be valid when 
ESI is turned off in Eq. (2). This term was neglected 
by Wolfenstein, which corresponds to the trivial solution 
of Eq. (21), assuming that the phase difference between 
the two CP-conserving amplitudes is zero for all decay 
channels. 

To be complete and detailing the notation of Ref. [2] 
by including the two-particle angular momentum states 
J, we show that the second term in Eq. (20), 




4 (sin7) E 


J 

OA 


A A' J 

i^OX ^a,a) ^OA 


( 22 ) 


also satisfies the CPT condition, namely, this quantity 
vanishes, which is easily verified by using Eqs. (17) as 


EAPf'^^ = 4 (sin 7 ) x 

A 


X 


E 

J 


XhXXJ (Box) ty x (^OA') 



(23) 


The vanishing of Eq. (23) is due to the symmetry of 
*A,A' = ^A'.A; and the fact that xxj = Xa'j, i-e., the 
strong interaction does not mix different CP eigenstates. 
Therefore, by taking into account Eqs. (21) and (23), one 
has that the CPT constraint 


E ATa = 4(sin7)Elm {Box}* A 

A A J 

+ EArf^' = o, 


(24) 


is fulfilled in leading order of the interaction. 


IV. RESONANT CHANNELS AND CPT 

In the case that the channel A contains also the for¬ 
mation of a resonance in the partonic process, namely. 


B —>■ 7r/5, the amplitudes Aqx and Box can be separated 
in the following two parts, = ^oAAffi+S/i Aqxr^ 

Bqx = Bqxnr+J^r Bqxrj where the subindex R and NR 
mean resonant and non resonant channels. Therefore, the 
decay amplitude in Eq. (19) is rewritten as 


^0XNR~ 


J 

p±*7 


^OXNR 


R 


A-lo “EE Aqxr + A 

i 

E ^OA'fi + Aox'NR 

R 

(E 


X',J 


'NR 


(25) 


The resonant source terms Aqxr and Bqxr should be 
interpreted as bare amplitudes, where at the resonance 
decay vertex, the two-hadron rescattering process is not 
yet included. The Breit-Wigner amplitudes for each term 
are identified according to 

(1 -h itAA)^0AK ^ al^FR^Pj(cos6) (26) 

and 

(1 + ^ti^)B^^^ ^ Pj{cos6l (27) 


where J is the spin of the resonance decaying to two spin 
zero particles and Pj {cos 9) is the Legendre polynomial 
and 6 is the helicity angle between the equally charge 
particles in the Gottfried-Jackson frame. We will give 
the representation of this angle for the B~^ —>■ tt+tt+tt” 
decay in Fig. 1 of the next section. 

After substituting (27) in (25), we get that 


Aio = E FiYPAcose) 

JR 

+ E (^OAATfl + Boxnr) 

J 

+ * E (Aox'nr + 'NR ) , (28) 

X'.J 


where the first and second terms in the right-hand side 
is the isobar model for the decay. The second term is the 
source term for the hnal state channel, and the third one 
includes the hadronic interaction among the two of the 
hadrons with angular momentum J. We should clarify 
that Eq. (28) includes the interaction in the resonance re¬ 
gion as the pair of hadrons has a probability to be formed 
directly from the partonic process. 
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The CP asymmetry from Eq. (28) can be cast in the following form 
APa = T{h^ X)-r(h^ A) 

= 4(sin7) 5]lm | Pj{cos9) + Bf ^ 


^OXNR 


Q\NR 


J J' 


OXNR 


+ * E E b^xFi^PAcos 6) + B, 

\> \ R 

- * E (E bl^yPi'yPj' (cos 0) + Bpy 


NR 


E ^oxFr'^ Pj' (cos 6) + Aqa 

/ \ R' 

^x',x ( E ^ox'FrN'Pj' (oos 6) + AA'n^ 

E ^oxFrx' PjA^s 0) + AAnr 


X' \ R' 


fF 

^X'.X 


(29) 


It has to be understood that the subindex A also in¬ 
cludes different kinematical regions of the three-body 
channel. Since we have introduced the Breit-Wigner am¬ 
plitudes in the decay amplitude, the CPT constraint has 
to be checked in the actual fit, i.e., if J2x = 0 is sat¬ 
isfied when one takes into account the integration over 
the phase-space besides the sum over all decay channels 
in the sum of A. Indeed, in our fitting procedure, we will 
keep only terms that, within our limited model, satisfy 
the CPT constraint. 


V. INTERFERING RESONANT AND NON 
RESONANT AMPLITUDES 


We present a simple example to explore the asymme¬ 
try formula (29), considering the resonant, non-resonant 
source terms and the contribution from the coupling be¬ 
tween two strongly interacting channels, namely tttt and 
KK. We use the vector and scalar resonances, the p(770) 
and /o(980) ones, for instance, interfering with a non res¬ 
onant amplitude and the term carrying the strong inter¬ 
action transition amplitude between the coupled chan¬ 
nels. This illustrates exactly the B^ —>■ de¬ 

cay case at low invariant 7r“'"7r“ mass. Also there are 
K^tt+ TT data [22, 23] previous to the CP asym¬ 
metry observation by LHCb collaboration. 

First, let us remind that, in a general way, the Breit- 
Wigner excitation curve for a resonance R reads 




1 

- s - imRTiiis) ’ 


(30) 


with mii being the resonance mass, and 


rfl(s) 


a-O''" 

(^-k) 


1/2 ’ 
gl/2 


(31) 


denoting the energy dependent relativistic width. For the 
pion mass, we adopted ttItt = 0.138 GeV, degenerated for 
the negative and positive charged particles. Here, we 
consider the resonance decay in the tttt channel. 


The real and imaginary parts of F]A^ (s) are given, 
respectively, by 


Re [F|W] = 

and 

Im [Fr 
The square modulus is 


mo — s 


oBWl 


{m% - sA + m%rR{sA' 
mRAnis) _ 

im\ - sA + m'Ar rAA ' 


(32) 


(33) 


\fR^\AA 


1 

(m| - sA + mj^TnisA ' 


(34) 


The amplitudes for B^ —>■ or 

^ K^Tr+ TT decays, taking into account the p(770) 
and /o(980) resonances interfering with a constant non 
resonant amplitude, can be written as 


= a^Ff'^kis) cos 0 + aiFf'^ 



+ %Ff'^k{s) COS0 + ( 35 ) 


where the kinematical factor k{s) = yl-^ is in¬ 

cluded in the amplitude of the p(770) vector resonance, 
to take into account the threshold behavior of the decay 
amplitude in a p-wave. The angle 0 is defined as the angle 
between the bachelor and the equally charged interacting 
particle. See this definition in the B'^ —>■ tt+tt+tt” decay 
illustrated in Fig. 1. Here, cos0 is associated to the spin 
1 of the p resonance and varies from —1 to -|-1 along the 
phase space. 



FIG. 1. —>■ 7r'''7r''"7r decay with tt'^ being the bachelor 

particle, (a): cos6 < 0 (0 > j). (b): cos0 > 0 (0 < 
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The form factor ^1 + j is included in the non 

resonant amplitude in order to parametrize the depen¬ 
dence on the square mass of the pair that is brought by 
the source terms, either for the tttt or KK systems. It 
comes from the partonic decay amplitude that produces 
the three-meson final state, in which the relative momen¬ 
tum between the pair of mesons is distributed among the 
quarks in the momentum loop within the microscopic 
amplitude, e. g., the tree and penguin diagrams, and 
probe the internal structure of the mesons involved in 
the initial and final states. In general, the non resonant 
amplitude should depend on the two kinematically inde¬ 
pendent Mandelstam variables, which for simplicity, we 
choose one of them. The suggestion is that the micro¬ 
scopic process gives raise to the momentum dependence 
in a power-law form reflecting the hard momentum struc¬ 
ture of the mesons involved in the decay [24] . For exam¬ 
ple, inspecting the tree diagram, one observes that the 
decay of the 6—quark in the rest frame, produces a fast 
light quark and a pion back to back. The light quark 
has to share its momentum with a slow antiquark, which 
is happens by the gluon exchange and is damped by a 
momentum power-law form. A similar reasoning can be 
applied to the penguin diagram. On the other hand, if 
factorization can be proved, the heavy to light general¬ 
ized transition form factors [25] like a heavy meson into 
iFTr/TTTr, for instance, will also play a role in CPV studies. 


An alternative to parametrize the decay amplitude, 

I 


convenient for Monte-Carlo simulations, is to write 
Eq. (35) as 


A'- 

*/-l( 




OA 


= cosOF -h 


(36) 


where and contain both the fixed weak and strong 
phases, with the Breit-Wigner functions introducing ad¬ 
ditional mass dependent strong phases as sketched above. 
The phase comes from the partonic amplitude pro¬ 
ducing the three-body final state, excluding the strong 
phase from the rescattering process. The relation be¬ 
tween the parameters is a^e*^± = , a£e*^± = 

al + and 

By comparing the first term in the r.h.s. of Eqs. (29) 
and (35), it is possible to make the following identifica¬ 
tions, namely, 


Aoxr = Wfc(s) cos 9 + a^Ff^, 

Boxr = cos 9 + blFf'^, 

nr 

Aoxnr = , I 

hnr 

Boxnr = ’ (37) 

^ + at 


and thus, rewrite Eq. (29) for the decay channels A = 
TTTTTT or Ktttt as 


ATa = 4(sin7)Im 


b^Ff^kis) cos 9 + b^Ff"^ + j j a^Ff^k{s) cos 9 + a^Ff^ + 


^ox 


I-f 


-I- 4(sin7) Re< ^ 


b^Ff'^kis) cos 9 + blFf'^ + 


hnr 

Cqa 


1 + 


j cos 9 + alFf'^ + - 


nr 

OX 


nr 

j.J=0 “oA' 

tx',x 




(38) 


where we used the identity Im(j 2 ;) = Re(z) and kept 
only one term in the scattering amplitude for J = 0 in 
Eq. (29). The subindex A are associated with a position 
in the phase-space of the B decay and A' is equal to ttKK 
or KKK, in the case of R —>■ tttttt and B — ^ Ktttt, re¬ 
spectively. The second term in the r.h.s. of Eq. (38) is 
the compound contribution to the CP asymmetry that 
includes the effect of the final state interaction in the 
non resonant channel. In this example, only the S'-wave 
amplitude KK ^ tttt channel was included and we have 
considered the interference of this term with the reso¬ 
nant ones. We use here the functional form of this term 
obtained in Ref. [2]. 


A. CP asymmetry formula 

The compound contribution to CPV, namely, the one 
carrying t'^Yx Eq. (38), has as input in our calculations 
the non-diagonal scattering amplitude depicted in the full 
S-matrix in the isoscalar and angular momentum zero 
TTTT KK coupled channels, which is written as 


zy/l — F’^s„„+&kk) 

ix/l-rf 


( 39 ) 
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with the inelasticity parameter 77 ( 5 ) and the tttt phase- 
shift given by [26], 


/ \ I k 2 — s , , 

?7(s) = 1 - ( ei^+ 62 —j -^-, (40) 


with 


fc 2 = 


y/s- 4to| 


K 


(41) 


respectively. In these expressions, we use rriK = 
0.494 GeV, M' = 1.5 GeV, Ms = 0.92 GeV, M/ = 
1.32 GeV, ei = 2.4, €2 = —5.5, and co = 1.3, according to 
the parametrization given in Ref. [26]. The off-diagonal 
term in the S-matrix gives the transition amplitude be¬ 
tween the TTTT and KK channels in the isoscalar channel, 
namely, = \/l — Tf k. y^l — 77 ^ ^ 

where we have made the approximation 5kk ~ be¬ 
tween 1 and 1.6 GeV. 


and 


'^TT-TT ('^) 


- COS 
2 


cot^[i5^^(s)] - 1 1 
cot^[(5^^(s)] + 1 J ’ 


with 


C 0 t(( 57 j- 7 I-) Cq 


{s - - s) \k 2 




q ’ 


The expression for the GP asymmetry in Eq. (38), 
should be expanded in a form where the unknown pa- 

(42) rameters are exposed in a simple manner to proceed with 
the fitting to the experimental data. By using the rela¬ 
tions ln\{ziz + Z 2 Z*) = Re( 2 :)Im( 2 i -|- Z 2 ) -I- Im(z)Re(zi — 
Z 2 ), Re(z*Z 2 ) = Re( 2 :i)Re( 2 : 2 ) + Im( 2 ;i)Im(z 2 ), and 
lm(ziZ 2 ) = Re(zi)Im( 22 ) — Ini(zi)Re(z 2 ): together with 

(43) Eqs. (33)-(34), one can finally write Eq. (38) as 


ArA= 


A 


{Scos[24^(s)] +B' sin[2(5^^(s)]} y/l - rf{s) 


+ |E®W(s)| 2 ^(s)cos 0 . 


(' + *) (‘ + 4 :) 


+ ClF^^(s)l^k^(s) cos^e 


V{ml - s) V'y/l - rf{s) { 7 nprp(s) cos[2(5^^(s)] - {ml - s) sin[2(5^^(s)]} 


1 -f 


-f 


1 -f 




fTOprp(s) £'y/T^^^ifJs) {{ml - s) C0s[2(5^,r(s)] + mpTp{s) sin[2(5,r7r(s)]} 


1 -f 


1 -f 




+ \Fjl^{s)\'^\Ff^{s)\'^k{s)cosex 

X {-^[("ip - s){m'j -s) + mpTp{s)mfTf{s)] 4- g[{ml - s)mfrf{s) - mpTp{s){m} - s)]} 

n{mj - s) 'H'yyi-'n^{s) |?7i/r/(s) cos[2(5^^(s)] - {mj - s) sin[2(5^^(s)]| 


BW/„m 2 




l + VT 1 +aI 


VmfTf{s) V'yJ^-Tf{s)< 

^^{mj - s) cos[2(5^^(s)] -f TO/r/(s) sin[2(5^^(s)]| | 


^ + ^ j 


+ Q|EfW(s)|2 


(44) 


where A is associated to or , and A' to 

tt^K'^K~ or K^K^K~. The masses of the resonances 
are mp = 0.775 GeV, and to/ = 0.975 GeV [ 8 , 27] . The 
widths used in Eq. (31) for R = p, f are T), = 0.150 GeV 
and Tj = 0.044 GeV. The parameters of the model can 
be written in terms of those in Eqs. (35) or (36) as follows, 

A = 4(sin7)Im = «a)^ - (a-A)^ (45) 

B = 4(sin 7 )Re - a^lh^ir ], (46) 

B' = -4(sin 7 )Im [a^j], b^l* + , (47) 


V = 4(sin7)Im [a^b'Sl* + 

= 2[aP^al\ cos(<5^ - J^a) - a'Lal\ cos(<5^ - T’^)], (49) 
P'=4(sin7)Im + <yK*] 

= 2lalal\, cos{6l - 6l\,) - cos{6^_ - 5rA/)](50) 

f = 4(sin7)Re - <xK*] 

= -2[alal\sin{Sl - Sl\) - a'Lal\sm{S^_ - (51) 

r = 4(sin7)Im [-a^y + 

= 2[alaTy sin(5^ - dl\,) + a>Lal\, sin((5^ - <5 ”(a0],(52) 


F = 4 (sin 7 )Im 


4K* + «n5n* 


0*^0 


= 2[a'{_a+ cos((5((. — (5:J() — a'Lai. cos((5(( — 5:(i)], (53) 


C = 4(sin7) Im [ag5g*] = {a^f - (a^()^ (48) 




























Q = 4(sin7)Re 


4 & 0 * “ “nfcn 




= 2[a^a{^_ sin(5^[. — 6^) — a'LaL sin(J^^ — St)], (54) 


TL = 4(sin7)Im 






= 2[aial\ cos{6i - 6^) - aLal\ cos(<5{ - r’,)], (55) 


the final state. In this case, the contribution from CPV 
comes only from the p(770) meson. In Fig. 2, we show 
the signature of the CPV from this term with the usual 
Breit-Wigner square modulus and with the same sign for 
both cases, namely, cos 9 < 0 and cos 9 > 0. This term 
locally violates the CPT constraint and will be set to zero 
in our fittings, by assuming C = 0. 


7^'= 4(sin7)Im Oq^ov + 

= cos(<5i - 6l\,) - aial\, cos(<5{ - (5rv)](56) 


P = 4(sin7)Re - a^lb^* 

= -2[aial\ sin(<5i - Sl\) - aLal\ s^sL - r(x)](57) 

P'=4(sin7)Im [-a^Xir + 0^,14* 

= sin{5i - 51\,) + sin(j{ - 5rv)],(58) 



and 


Q = 4(sin7)Im 


fuf* 

^0% 


(59) 


In the previous work of Ref. [2], we used a differ¬ 
ent form for the KK —>■ tttt amplitude, writing it as 
IRTaI cos(^a + Sy + 4'a), where Kx = B^^Aoy - BoxA^y 
and <I)a = —i ln(RrA/|R^A|)- This form is analogous of that 
used in Eq. (44), in terms of B and B'. The phase 4 >a 
was chosen to be zero that is the same as to take B' = 0, 
assumed hereafter. 

As a remark, we call the reader attention that for the 
compound contribution, we have discarded three-body 
rescattering effects at the two-loop level since they are 
small compared to the first two-body collision contribu¬ 
tion, as suggested by the three-body model calculations 
for the decay, see Refs. [28-30]. 


VI. ANALYSIS OF CP ASYMMETRY TERMS 

In order to be able to compare our model di¬ 
rectly with the recent experimental data regarding the 
—>■ 7r=*=7r+7r“ and B^ —>■ TT decays presented 

in Ref. [9], we had to eliminate one variable in Eq. (44), 
namely, the square mass of the pair presenting the bache¬ 
lor particle. We perform this calculation in Appendix A. 
Here in this section, we show the behavior of some terms 
of the integrated Eq. (44) . 


A. Direct CPV term containing the C parameter 

As it was mentioned before, the term of the integrated 
Eq. (44) containing the C parameter, is equivalent to the 
direct CPV induced by the interference between the tree 
and penguin amplitudes with one meson resonance in 


FIG. 2. Term containing the C parameter of the integrated 
Eq. (44), as a function of m^+^- = ^/s for (a) cosO < 0 and 
(b) cos 9 > 0. In these plots we set C = 1. 


B. DCPV term containing the D parameter 

The Dalitz CPV induced by the interference term of 
the integrated Eq. (44), namely, that containing the V 
parameter, is directly related with the real part of the 
Breit-Wigner function. This function presents a clear sig¬ 
nature in the mass spectrum, namely, a zero and a change 
in the sign of the CP asymmetry at the central value of 
the p(770) mass. Another peculiarity associated to vec¬ 
tor mesons resonances is one more sign change when the 
COS0 cross the zero around the middle of the Dalitz plot. 
These two facts can be observed in Fig. 3. Notice that 
this term sums to zero and thus locally satisfies the CPT 
constraint. 


C. DCPV term containing the £ parameter 


The term containing the £ parameter in the integrated 
Eq. (44), associated to the Dalitz CPV induced by the 
interference between the p(770) and the non resonant am¬ 
plitude, is directly related with the imaginary part of the 
Breit-Wigner function of p(770). The shape is similar 
to the those presented in Eig. 2, with the clear difference 
that the proportionality with cos 9 changes the sign of the 
CP asymmetry when cos 9 pass through zero in the mid¬ 
dle of the Dalitz plot. Fig. 4 shows these features. Notice 
that this term sums to zero and thus locally satisfies the 
CPT constraint. 
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FIG. 3. Term containing the T> parameter of the integrated 
Eq. (44), as a function of for (a) cos 9 < 0 and 

(b) cos 6 > 0. Here we have used A.\ = 1.705 GeV and 11 = 1. 



FIG. 4. Term containing the £ parameter of the integrated 
Eq. (44), as a function of = ^/s for (a) cos 9 < 0 and 

(b) cos 9 > 0. Here we have used = 1.705 GeV and £ = 1. 


D. DCPV term containing the J" parameter 

The projection of the difference between B~ and B'^ 
events presents a clear signature in the mass spectrum 
with a zero close to the mass of the vector meson p, and 
another one close to the mass of the scalar /o(980). There 
is also a change of sign in the CP asymmetry associated 
to the cos 9 passing through zero around the middle of 
the Dalitz plot. These features can be observed in Fig. 5, 
where we display the term containing the T parameter 
of the integrated Eq. (44). This term locally satisfies the 
CPT constraint. 


E. DCPV term containing the Q parameter 

The shape of this term presents the characteristic 
peaks of the two interfering resonances. Due to the direct 
proportionality to cos 6, there is a change of sign of the 
CP asymmetry when cos 6 pass through zero in the mid¬ 
dle of the Dalitz plot. Fig. 6 shows the features of this 



FIG. 5. Term containing the F parameter of the integrated 
Eq. (44), as a function of m^+^- = ^ys for (a) cos 9 < 0 and 
(b) cos 9 > 0. Here we set = 1. 

term. This term is consistent with the CPT constraint 
locally. 



FIG. 6. Term containing the G parameter of the integrated 
Eq. (44), as a function of m^+^~ = y/s for (a) cos 9 < 0 and 
(b) cos 9 > 0. Here we set G = \. 


VII. CPV IN 7r±7r+7r- AND ->■ tt^A+A- 

DECAYS 

We determine the relative contributions related to the 
integrated Eq. (44) by fitting two CP asymmetries dis¬ 
tributions presented in Ref. [9]; one for the tt+tt” mass 
distribution coming from the difference between B~ and 
in —>■ TT^TT+TT decay, in particular the distri¬ 

bution with cos 9 > 0. From this fitting we found the 
parameters and plot the CP asymmetry for cos0 < 0. 
The other fit was performed for the phase-space integra¬ 
tion in the region cos 9 < 0 of the tt+tt" mass distribu¬ 
tions in the B^ —i Ar=*=7r+7r“ decay. From these two hts 
we get the value of the tt+tt” —>■ K^K~ parameter {B) 
for each decay channel and plot the asymmetry for the 
K'^K~ distribution in the channels B^ —>■ •k^K^K~ 
and ^ K^K+K-. 
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To perform both fits and the other plots, we use the 
7r+7r“ —>■ K~^K~ amplitude, the non resonant compo¬ 
nent, and the p and /o(980) resonances. We put these 
amplitudes within the isobar model through a coherent 
sum of them. The biggest source of uncertainty from 
these fits is the error in the tttt phase shift and inelastic¬ 
ity parameter. For simplicity we use the central value of 
Ref. [26], although possible variations within the quoted 
errors in the parameters could give better results. De¬ 
spite our belief that this model is able to represent much 
of this data, we are aware that it can not explain all the 
rich CP violation structure observed in these decays. The 
inclusion of other contributions and the symmetrization 
of the —?► 7r=*=7r+7r“ decay amplitude would be nec¬ 

essary to understand in more details the phase space of 
these decays and get better agreement in all CP asym¬ 
metry regions presented in Ref. [9]. 

To perform the fits we use for the /o(980) resonance, 
ruf = 0.975 GeV and width P^ = 0.044 GeV, getting 
from E791 experiment [27]. For the two A parameters, 
we use Aa = A^^r = 3.0 GeV and Ay = Akk = 4.0 GeV. 
It is important to say that we changed these values by a 
factor two without finding an appreciable change in the 
fitting. 


A. TTTrTT channel 

We started by performing the fit of the integrated 
Eq. (44) to the cos 0 > 0 asymmetry distribution of the 
low 7r’''7r“ mass projection for the —>■ decay. 

Due the presence of identical particles in hnal state, the 
interference between symmetrical terms by the exchange 
of the two identical pions must disturb the CP violation 
pattern mostly for cos 9 < 0. We study the regions on 
the Dalitz plot where this interference can be minimum 
and observe that it corresponds to the high tt+tt” region 
on the cos 9 > 0 distribution as can be seen in Fig. 7. 



FIG. 7. as a function of cos 6 for some particular 

values of m^+^-. 

This figure provides the mass of the changeless pair of 
pions formed with the bachelor pion. If cos 9 < 0 the 
mass of this pair stays below 3 GeV, and for < 1 

GeV, it can be even below 1 GeV, and therefore mak¬ 


ing relevant the interference with the resonances. For 
cos 9 > 0, my+^- > 3 GeV, minimizing interference 
effects from the Bose symmetrization of the decay am¬ 
plitude. We intend to perform this study in the future. 
For the moment, we do not consider the symmetrization 
of the decay amplitude, and therefore we do not use the 
data for the GP asymmetry in the B^ —>■ chan¬ 

nel for cos 0 < 0 to fit the parameters. 

The best fit was obtained by using four parameters 
associated to the tt+tt” —>• K'^K~ amplitude [B), the 
real and imaginary parts of the interference between the 
p and the non resonant partonic amplitudes {B and f), 
and finally the imaginary part of the interference between 
the p and the /o(980) resonances [Q). The result is pre¬ 
sented in Fig. 8a, where we compared our model with the 
experimental data extracted from Fig. 4c of Ref. [9]. The 
remaining parameters, namely, those that locally violate 
the GPT constraint, and those negligible in the fitting 
procedure, were set to zero. 



FIG. 8. (Color online) CP asymmetry of the —>■ 7r^7r’'"7r“ 

decay, integrated Eq. (44), compared with the experimental 
values (blue points) taken from Ref. [9]. Results for cos 9 > 0 
for (a) total and (b) individual contributions. 

The individual contributions to the GP asymmetry are 
shown in Fig. 8b. As expected, the p meson contribution 
represented by the amplitudes containing the V and £ 
parameters, is large mainly for the real part of the BW, 
and the presence of /o(980) is seen only by one of the 
interfering terms with the p, namely, that presenting the 
Q parameter. All these terms are locally GPT invariant, 
namely by integration in cos 9 their contribution to the 
GP asymmetry vanishes. Between 1 and 1.6 GeV, the 
contribution of the tt+tt” —^ K^K~ amplitude, namely, 
the B term, shows its importance and as expected [2], 
dominates the asymmetry in this region. We remind 
that this part of the asymmetry does not vanishes upon 
integration in cos 9 and cancels the asymmetry in the 
B^ —>■ tt^K^K~ decay channel. 

Looking at the behavior of the experimental cos 0 < 0 
distribution in Fig. 9a, we can see that much of the fea- 
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tures observed for the distribution of events for cos 0 > 0 
are present in this asymmetry. In fact, we plotted to¬ 
gether with the experimental points, the amplitudes com¬ 
puted with the integrated Eq. (44) using the fit param¬ 
eters obtained from the cos 9 > 0 distribution. The 
plot in Fig. 9a shows a clear departure from the exper¬ 
imental data when at the starting of the p mass reso¬ 
nance and also at the beginning of the contribution of 
the 7r+7r“ —>■ K~^K~ amplitude to the asymmetry. As¬ 
suming that these differences are due to the symmetriza- 
tion, this result suggests an interference between p and 
7r+7r“ —^ K~^K~ amplitude in the crossing channels, 
which is corroborated by Fig. 7, where for the mass region 
above 1 GeV, the mass of the tt'+tt” pair in the crossing 
channel for cos 6 < —0.75 can be even below 1 GeV. In 
addition, we show in Fig. 9b the individual contributions 
to this GP asymmetry distribution. 



FIG. 9. (Color online) CP asymmetry of the —>■ 7r^7r"'"7r“ 

decay, integrated Eq. (44), compared with the experimental 
values (blue points) taken from Fig. 4d of Ref. [9]. Results 
for cos 0 < 0 for (a) total and (b) individual contributions. 

We still remark to the reader that for the —>■ 

7r^7r+7r“ decay, our definition of the angle 9 is opposite 
from that presented in Ref. [9]. Thus, we had to compare 
our results of cos 0 > 0 (cos 0 < 0) with those of cos 9 < 0 
(cos 9 > 0) of Ref. [9]. 


B. -kKK channel 


The expression for the CP asymmetry in the coupled 
^ channel, derived from the general for¬ 

mula in Eq. (29), applied to this specific decay and inte¬ 
grated in cos 9, is given by 


Ar(s) 


_2A_ 

a'{s)y/s -(l -h 

2By^l - p^js) cos[2J^^(s)] 
a'(s)\/s-4m^ (l + 7^) (l + 7 ^) 


where the kinematical factors, namely. 


a'(s) = 


(s — 


^ B ^ ^ _ ^2 


4s 


1/2 ’ 


(61) 


and the kaon momentum in the rest frame of the KK 
subsystem, s — 4:171^, are now taken for the ttKK sys¬ 
tem. Furthermore, the integrated decay width in cos 9 
from Eq. (60), becomes exactly opposite in sign to the 
decay width in the tttttt channel above the KK thresh¬ 
old, c.f. the first two terms of the integrated Eq. (44). 

The expression in Eq. (60) is obtained by adding the 
7r+7r“ —>■ K'^K~ contributions detailed in Figs. 8 and 
9, since A = 0. The plot shown in Fig. 10 was done 
with the parameters fitted by the CP asymmetry data 
in the B^ —>■ 7r=*=7r+7r“ decay. The CP violation dis¬ 
tribution obtained for the K^K~ invariant mass from 
B^ —>• tt^K^K~ decay through Eq. (60), has opposite 
sign with respect to the correspondent term in the cou¬ 
pled 7r^7r''"7r“ channel. 



FIG. 10. (Golor online) CP asymmetry of the 
—>■ tt^K^ K~ decay, Eq. (60), compared with ex¬ 
perimental data (blue points) taken from Fig. 7b of Ref. [9]. 


For this decay, the LHCb experiment presented the 
sum of events with cos 0 < 0 and cos 9 > 0. The compar¬ 
ison of the data with our model is also provided in Fig. 10. 
There is a clear agreement about the shape of the model 
distribution with the experimental data and also a rea¬ 
sonable amount of the number of events related with this 
kind of CP asymmetry. It is also worth to mention that 
CPT is conserved if we sum all CP asymmetry contri¬ 
butions obtained with our approach to B^ —>■ 7r^7r''"7r“ 
and B^ -y 'k^K^K~ decays, in a region of tt+tt” and 
K'^K~ invariant mass below 1.6 GeV. 


VIII. CPV IN K^'k+tt- and 
B± K^K+K- DECAYS 

The number of events observed in these two decay 
channels are about one order of magnitude larger than 
in the and Tr^K'^K~ data [9]. This allows a 

better fit to the difference between B~ and B'^ events. 
In the fitting procedure we use only the cos 0 < 0 distri¬ 
bution, because the experimental results [9] for cos 0 > 0 
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present a new feature in both decays in the region stud¬ 
ied in this work. Our simple model does not account for 
this new behavior. We only have a guess, that we already 
mentioned in our previous paper [2] , related with the pos¬ 
sible presence of re-scattering coming from double charm 
decays. One should note that three light-pseudoscalar 
mesons can, in principle, couple via strong interaction 
with channels like DDK. It seems reasonable to expect 
that DDK —>■ KKK or KinT can contribute to the CP 
asymmetry in regions of large two-body invariant mass 
above the DD threshold, that is far from the KK thresh¬ 
old and above 1.6 GeV, outside the region discussed in 
this work but excluded from the cos 0 < 0 distribution. 


comes from the real part of the interference between p 
and the /o(980) resonances. The imaginary part of this 
amplitude has a contribution, however, much less impor¬ 
tant than the real one. The 7r“*"7r“ and K~^K~ amplitude 
plays an important role in this fit above 1 GeV and also 
associated with the interference with the p resonance. 

Here, we also used Aa = 3.0 GeV and A^/ = 4.0 GeV. 
The other parameters were found by using the method 
in order to fit the experimental data for cos 0 < 0 distri¬ 
bution from Ref. [9]. 

B. KKK channel 


A. Ktitt channel 


Differently to the —?> 7r=*=7r+7r“ decay where the 

p amplitudes were dominant, here the /o(980) has the 
largest contribution to the ^ AT^tt+tt" decay, as it 
is expected by recalling the BSS mechanism applied to 
these decays, which builds in principle the correspond¬ 
ing amplitudes associated with the source terms or par- 
tonic amplitudes. One possible good fit for the integrated 
Eq. (44) is obtained by using five parameters associated 
to the TT+TT” —>■ K^K~ amplitude (H), the real parts of 
the interference between p resonance with the non res¬ 
onant partonic (D), and tt+tt” —^ K^K~ amplitudes 
{D'), and finally the real and imaginary parts of the in¬ 
terference between the p and the /o(980) resonances, {T 
and Q). As in the B^ —>■ case, the parameters 

that locally violate the CPT constraint, and those negli¬ 
gible in the fitting procedure were set to zero. Fig. 11a 
shown the best fit we get to the cos 6* < 0 distribution. 




The GP asymmetry in the B~ and B^ event distri¬ 
butions in the B^ —>■ K^K'^K~ channel is compared 
to our model with the sum of the events of cos 9 < 0 
and cos9 > 0 as we did for B^ —>■ Tr^K'^K~ decay 
channel and the data provided by the LHGb in this last 
case. The functional form of the asymmetry for the 
B^ —>■ K^K~^K~ decay is the same as in Eq. (60). 
The only change is due to the kinematical factor a'{s) in 
Eq. (61), that now presents the replacement —>■ 

We plot in Fig. 12 the asymmetry of the B^ —>• 
K^K~^K~ decay, computed with the parameter B ob¬ 
tained in our previous fitting of the B^ —>■ de¬ 

cay. Our model is compared to the data obtained by 
summing the experimental results of B^ —>■ K^K'^K~ 
distributions for both cos 9 positive and negative regions. 



FIG. 12. CP asymmetry of the —>■ K~ decay com¬ 

pared with experimental values (blue points) taken from the 
sum of Figs. 6c and 6d of Ref. [9]. 


IX. CONCLUDING REMARKS 


FIG. 11. (Color online) CP asymmetry of the —>• 

A'^7r’''7r~ decay, integrated Eq. (44), compared with the ex¬ 
perimental values (blue points) taken from Fig. 5c of Ref. [9]. 
Results for cos6 < 0 for (a) total and (b) individual contri¬ 
butions. 

In Fig. 11b, we shown the decomposition of each com¬ 
ponent of the fit. Clearly the dominant contribution 


A. Beyond the model 

We can go beyond the asymmetry formula in Eq. (38), 
where we have not included the interference terms be¬ 
tween the resonances and the scattering terms in the elas¬ 
tic channel, we remind that in the BW we included the 
elastic amplitude, but it is also important to consider it 
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outside the resonance region, and returning to (29) as the 
starting point. We have neglected the Bose symmetriza- 
tion of the amplitude by the exchange of the identical 
pions and the contribution of the double charm produc¬ 
tion and the its coupling to the KK channel. This also 
should be addressed in the future. 

The isoscalar and zero angular momentum tttt 
originated by contribution of the elastic scattering 
terms to the CP asymmetry for —>■ 7r^7r“*"7r“ and 

—>■ K^Tr'^7T~ decays written before, can be substi¬ 
tuted the inelasticity and phase-shifts given by tttt —>■ 
KK S-matrix provided by Eq. (39) as = 

i (l — The inclusion of the S'—wave tttt scat¬ 

tering should improve the fitting below the p, where it is 
clearly seen in Fig. 9, the effect of a missing interference 
term. In this particular region the /o(600), namely the 
cr, will be taken into account by and we expect to 

improve the fit in the low mass region. 

The interference between tttt —^ ifiC and the p reso¬ 
nance was used to fit the Ktttt channel, as we can see in 
Fig. 11. Other possibility in this fit would be consider 
the large presence of /o(980), by using its interference 
with TTTT —>■ KK. This component do not satisfy CPT 
locally, but is consistent with Eq. (29) with the coupled 
channel B^ —>■ K^K^K~. In this work, we choose to 
use only locally CPT invariant components for the inter¬ 
ference with resonances. 


B. Conclusions 

We show how CP asymmetry in charmless three-body 
B^ decays is constrained by CPT in the presence of res¬ 
onances and final state interaction in leading order. The 
CPT constraint AP^ = 0, should include in the sum 
over the final states, all the kinematical allowed phase 
space and possible decay channels. In leading order, the 
scattering matrix corresponding to the decay chan¬ 
nel where a resonance is formed, e.g., the p in p-wave tttt 
elastic amplitude, accounts for the dressing of the decay 
vertex p — >■ tttt. The physical decay amplitudes associ¬ 
ated with the partonic source, that give the p meson and 
in its sequential decay to two pions, is dressed by the 
TTTT re-scattering process in the resonant channel. The 
non-resonant channels receive contribution from the re¬ 
scattering in leading order. The decay of the resonance 
and re-scattering by the bachelor meson, should appear 
as a next-to-leading order contribution, to be consider in 
future works, where the general formulas for the three- 
body decay amplitude based on the CPT constraint that 
includes the formation of resonances, re-scattering and 
coupling to inelastic channels should be addressed. In 
our limited two-body description, we explored the rich 
patterns of the CP asymmetry formula, namely interfer¬ 
ence between non-resonant and resonant terms, as well 
as the coupling between the tttt and KK channels due to 
the strong interaction. 

We analysed in our study two examples of CP viola¬ 


tion in charmless B^ three-body decays in the presence of 
resonances, excluding the possibility of channel coupling. 
The first case discussed is the p(770) meson resonance 
plus a non resonant channel, and the second one is the 
case of two interfering resonances, namely, p(770) and 
/o(980). These examples discussed within two amplitude 
models are suitable for the study of B^ —>■ 7r=*^7r+7r“ and 
B^ —^ i^^7r+7r“ for masses below the KK thresh¬ 

old, as provided by the new data LHCb data [9] , and we 
selected terms that are locally constrained by CPT. 

The CP asymmetry produced by the p(770) meson res¬ 
onance in the simple two amplitude model with a non¬ 
resonant channel, which are brought by the real and 
imaginary parts of the interference term, presents a clear 
signature in the low mass tt+tt"'' spectrum for both de¬ 
cays B^ —^ 7r^7r“''7r“ and B^ —^ K^tt~^tt~, as indeed 
the new data suggests [9]. One term shows a very clear 
signature: a zero at the p(770) mass in both regions of 
cosd < 0 and cos0 < 0 for B^ —>■ tt^tt+tt” with a change 
of the sign around cos 0 = 0. 

The two amplitude model constituted by the reso¬ 
nances p(770) and /o(980), shows two zeros in the CP 
asymmetry from the real part of the Breit-Wigners: 
one close to the p(770) mass and another one close to 
the /o(980) mass. However, the term which mixes the 
real and imaginary parts of the Breit-Wigners does not 
present a clear signature of CP violation around the 
p(770) meson mass, while a zero is present at the /o(980) 
mass, because the ratio Tpirip/Tf rrif is considerably 
larger than one. In both examples, the zero’s of the asym¬ 
metry at the resonance mass are associated only with the 
Dalitz CP violation from the interference of two ampli¬ 
tudes and the real part of the Breit-Wigner, the contribu¬ 
tions associated with the imaginary part present a behav¬ 
ior with a maximum at the resonances masses. The re¬ 
cent experimental results [9] are strongly indicating that 
the CP asymmetry appearing in the low mass re¬ 

gion comes from the real part of the Breit-Wigners in 
the interference between two amplitudes involving the p 
meson resonance. 

In addition, the CP asymmetries in whole kinemati¬ 
cal region below rriT^+T^- < 1.8 GeV in B^ 7r=*=7r+7r“ 
and B^ —>■ K^tt~^ TT decays demand the coupling to the 
7r^7r“*"7r“ by the strong interaction, as has already been 
discussed in Ref. [2]. To account for that, we introduced 
in the fitting of the CPV, the contribution of the coupled 
TTTT to KK channels in isoscalar and zero angular states. 
This was relevant for the fitting above the KK thresh¬ 
old. The independent confirmation of the CP asymmetry 
flowing between different channels coupled by the strong 
interaction was checked by the fair reproduction of the 
CPV in B^ —>■ Tr^K^K~ and B^ —>■ K^K^K~ distri¬ 
butions. The computation of the asymmetry was per¬ 
formed as dictated by the CPT constraint applied to 
coupled channels. There was no new parameters beyond 
those given in the fitting of the CP asymmetry data in 
the B^ —>■ 7r^7r+7r“ and B^ —>■ K^Tr~^Tr~ decays. 

We discussed improvements to CP asymmetry formula 
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beyond the one used to fit the experimental data on 
^ TT^TT+TT” decay, where the interference between 
the resonances and elastic scattering terms are missing. 
The inclusion of the interference between the S'—wave tttt 
scattering amplitude and the p(770) meson Breit-Wigner 
should improve the fit below the rup. In this particular 
region, the /o(600), namely the tr, will be taken into ac¬ 
count through the elastic tttt scattering, which it is ex¬ 
pected to improve the fit. Still the Bose symmetrization 
of the amplitude by the exchange of the identical pions 
has to be considered, which should be addressed in the 
future. 

Finally, the results of our fittings suggest an impor¬ 
tant consequence to CP violation: the strong phases can 
come from the hadronic interaction minimizing the effect 
of the imaginary part of penguin contributions to both 
decays studied here, as the recent experimental results 
are indicating. 


and 


b{s) = -- 


Mo -I- 3m^ — s 


2(s- 4m2)i/2 


(Mr— s)2 r) 

^ B ^ ■ - ml 


-, 1/2 


As 


, (A3) 


with Mb = 5.279 GeV (the detailed derivation of this 
equation is relegated to the Appendix B). Therefore, the 
asymmetry presented in Eq. (44) is actually a function of 
s and 77132, due to cos 9. In order to eliminate m ‘^21 
tegrate ArA(s, m‘ 12 ) in this variable, leading to the asym¬ 
metry Ar(s) = / ArA(s, m§2)c^7T7|2- We still separate 
the integration in two regions, namely, that of cos9 < 0, 
leading to 


Ar(s)(™«e<o) 



ArA(s,m^2)c^™32> 


(A4) 


Appendix A: Angle integrated asymmetry 

By defining in Fig. I, 7r+ = I, 7r“ = 2, and = 
3 (bachelor), one has that m^+^- = 77112 = a/s, and 
m^i+^- = 77132. The kinematics of the decay imposes 
that cos 0 is a function of s and 77132 , namely, 


and that of cos0 > 0, leading to 

Ar(s)(™®®>°) = / ArA(s, 77132)d77l32. (A5) 

J —hj a 

The quantities 


COS 0 = 0(5)77132 -f 6 (s), 


(Al) 


where, 


a(s) = 


(5 — 4771^)^/^ 


[Ml-ral-sY _ 2 

- As - ilhL 


1/2 


, (A2) 


(n7y-i(5) =(m2 = L^(A6) 

a(sj a(s) 

correspond, respectively, to cos 6 = —I and cos0 = -l-I, 
according to Eq. (Al). The quantity —b/a is related to 
COS0 = 0. 


The integration in Eqs. (A4) and (A5) generates the following asymmetry. 


Ar(s)= 


A 


Bcos[2S^^is)Wl-r,^s) , {s)\^k^{s) 


i{s)^s-4ml(l + ^y a{s)^s - 4ml (l + ^) (l + 


3a(s) 


-h 


|F/'^(s)|2fc(s) j T>(77i 2 - s) 27 'a/ 1 “ {"iprp(s) cos[2^^^(s)] - (m^ - s) sin[2^^^(s)]} 


2^a(s) 


1 + 


\/s - 477 i2 (l + 

SmpTpjs) g'A/l-7?2(5){ (?7ig - s) cos[2(5 TTTT (s)] -I- ?7ipFp(s) sin[25^^(s)]} 
1 + * ^ 




|ABW(,)|2|ABW(,)|2fe(^) 

2^a(s) 


X {F[{ml - s)(jn) - s) -h 771pFp(s)771/F/(s)] -t- a[(771p - s)771/r/(s) - 771pFp(s)(771^ - s)]} 
|F/W(s)|2 J -^(^2 _ g) 'H'A/l-f? 2 (s) |771/F/(s) cos[2^^^(s)] - {m} - s) sin[2(5^^(s)]| 


a(s) 


I-f 


\/s - 477 i2 (1 + 


rm/Ffis) {("»/ - s) cos[2(5^^(s)] -h mfTf(s) sin[2(5^^(s)]| 1 Q|E/W(s)|2 


1 + 


a/s - 4?7i2 (i -H 


i(s) 


(A7) 
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where ^ = — 1 for cos 0 < 0 , and ^ = 1 for cos 0 > 0 . 
Here, the rest frame momentum of the pion, namely, 
(s — is included in the terms presenting in¬ 

terference with the TTTT KK amplitude, in order to 
introduce the kinematical factor in the scattering ampli¬ 
tude. Furthermore, this factor is consistent with the CPT 
constraint, which indicates that the decay widths for the 
coupled TTTTTT and ttKK channels above the KK thresh¬ 
old is such that AF^ttii- = — AF^ricic, as discussed at the 
end of Sec. VIIB. Finally, we also remind that all free 
parameters are compatible with the energy quantities of 
the model, that are in GeV units. 


Appendix B: Derivation of Eq. (Al) 


In a general way, one has that 

r)U T)f 0/| 1/| 2/| 3/ 

P P^J. =P Po+P Pi +P P 2 +P Ps 

= EE' - Pa;p'^ - Pyp'y - p^p'^ 

= EE'-p-p'= pop'q-p-p'. (B1) 

If then Eq. (Bl) becomes 

= E"^ — p"^ = + P^ — P^ = ■ (B 2 ) 

In the —>■ tt+tt+tt" decay of Fig. 1, we define 7 r+ = 

1, 7 r“ = 2, and ■n''^ = 3 (bachelor). Therefore, = 

TOi 2 = \/s, rriT^i+T^- = 17132 , and rrii = m 2 = m 3 = m,r, 
since we have degenerated the pion mass. If we assume 
that 


-pl^ 

-^32 


= p^ + p^, 


(B3) 


then, these definitions lead to 

mh = P^2 P32^ = {P3^ + P^){P3f. + P2^) 

_ pM p _i_ pM p I pM p _i_ pM p 

— J 3 -^3 M -^3 -^2 + -^2 -^3 /i “T -^2 -^2 

= ml + P 3 P 2 -P3-P2+ P 2 P 3 -P2-P3+ml 
= 2ml + 2 pIpI - 2p3 ■ p2 
= 2ml + 2 pIpI + 2\p 3\\p 2 1 cos0. (B4) 


Thus, cos 9 is written as 


cosd 


"^32 - - 2pIpI 

2|p3||p2| 


(B5) 


In order to write cosd in terms of TO 32 and m ^2 = 
it is necessary to find P 2 , \P 2 \, P 3 , and jpaj. The first two 
quantities can be extracted as follows. Notice that, 

TO?2 = S = P^2 ^12 M = [Pi + P2){Pl,^ + 

_ p/^ p _i_ pM p _i_ p/^ p I pM p 

— 1 1 /.A "T J 1 j 2 2 -^1 M -^2 -^2 /.A 

= ml+ p\p 2 - Pi ■ P 2 + P 2 P 1 -P 2 -Pi+ml 
= 2ml + 2 pIpI - 2pi ■ p2 
= 2to 2 + 2p?p^-t 2 |pi||p2|- (B 6 ) 


Furthermore, by the figure, we see that |pi| = |p 2 | = |Pir|- 
However, since p\ = yjm\ -|- |pi p and P 2 = i/mf -I- |pip, 
we conclude that p\= P 2 =p^ = -I- IpVp- Then, 

s = 2ml + 2{ml + Ip^l"^) + 2|p^p = 4m^ -b 4|p^|^ (B7) 
with 

and 

p° = \/ml + IpttP = (B9) 

In order to find the s dependence of pg and jpaj, we 
proceed to write 

P^ = Pf -b P^ + P^ = Pf 2 -b P^, (BIO) 

and make 

Ml=P^PBf. = {Pr2 + P3 ^KPi2i. + P3f.) 

= Pf2Pl2f. + Pi^ 2P3>. + P^Pl2>. + P^P3^ 

= s+ P? 2 P 3 - P 12 ■ P 3 + pIpi 2 - P3 ■ P 12 + ml. 

(Bll) 

In the reference frame in which pi 2 = pi -b P 2 = 0, we 
have 


Ml = s-b2p?2P3 +m^. 
Furthermore, in this reference frame. 


\/s+ IP 12 P = y/s, leading to 

MI = S -b 2 '\/sp 3 -b ' 
From this equation we get, 


- 
P3 — 


Ml-ml-. 


2 yi 


(B12) 

P 12 


(B13) 

(B14) 


Finally, from pg = m| -b |p 3 p, we have 

IpsI = \J{plf - ml 


{Ml -ml- sY 
45 


(B15) 


Inserting Eqs. (B 8 ), (B9), (B14) and (BI5) in Eq. (B5) 
we get. 


COS0 = 


™32 - - 


Mr— mi.—s 




1/2 


(M^—m^ — s )2 


4s 


1/2 


2777.32 — M% — 3ml + s 


2{s — 4 m 2 )^/^ 


(Ml-ml-sY 


-, 1/2 


4s 


. (BI 6 ) 


This same equation is obtained if we now treat the 
B~ —>■ 7r“7r+7r“ decay. In this case, the B meson and the 
pions of Fig. 1 have its charges changed, and one defines 
7r“ = I, 7r+ = 2 and 'k'~ = 3 (bachelor). Therefore, one 
has m^r+Tr- = m 2 i = ^/s, and m.^+.^,- = 77123 . For this 
case, cosd will be written exactly as in Eq. (B16), but 
for 77732 —>■ 77123 . 
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